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The two main advantages of space-based observation of extreme energy (> 5 x 10^® eV) cosmic rays 
(EECRs) over ground based observatories are the increased field of view and the full-sky coverage 
with nearly uniform systematics across the entire sky. The former guarantees increased statistics, 
whereas the latter enables a clean partitioning of the sky into spherical harmonics. The discovery of 
anisotropies would help to identify the long sought origin of EECRs. We begin an investigation of 
the reach of a full-sky space-based experiment such as EUSO to detect anisotropies in the extreme- 
energy cosmic-ray sky compared to ground based partial-sky experiments such as the Pierre Auger 
Observatory and Telescope Array. The technique is explained here, and simulations for a Universe 
with just two nonzero multipoles, monopole plus either dipole or quadrupole, are presented. These 
simulations quantify the advantages of space-based, all-sky coverage. 


I. INTRODUCTION 

A. Previous Anisotropy Searches 

No all-sky observatory has yet flown. Consequently, 
the first full-sky large anisotropy search was based on 
combined northern and southern hemisphere ground- 
based data. The respective data came from the 
SUGAR [T] and AGASA [5] experiments, taken over 
a 10 yr period. Nearly uniform exposure to the en¬ 
tire sky resulted. No significant deviation from isotropy 
was seen by these experiments, even at energies beyond 
4 X 10^® eV [5]. More recently, the Pierre Auger Gol- 
laboration carried out various searches for large scale 
anisotropies in the distribution of arrival directions of 
cosmic rays above 10^® eV (an EeV) [H [S]. At the en¬ 
ergies exceeding 6 x 10^® eV, early hints for a dipole 
anisotropy existed, but these hints have grown increas¬ 
ingly weaker in statistical strength [6j. The latest Auger 
study was performed as a function of both declination 
and right ascension (RA) in several energy ranges above 
10^® eV. Their results were reported in terms of dipole 
and quadrupole amplitudes. Assuming that any cosmic 
ray anisotropy is dominated by dipole and quadrupole 
moments in this energy range, the Pierre Auger Collab¬ 
oration derived upper limits on their amplitudes. Such 
upper limits challenge an origin of cosmic rays above 10^® 
eV from long lived galactic sources densely distributed in 
the galactic disk [7]. In the E > 8 EeV bin, they did re¬ 
port a dipolar signal with a p-value of 6.4 x 10“® (not 
including a “look elsewhere” penalty factor) [8]. Their 
cutoff of ~ 8 EeV is above the galactic to extragalactic 
transition energy of ^ 1 EeV, but still below the GZK 
cutoff energy of ~ 55 EeV. Also, Telescope Array (TA), 
the largest cosmic ray experiment in the northern hemi¬ 
sphere, has reported a weak anisotropy signal above their 
highest energy cut of 57 EeV [H]. 
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B. Extreme Universe Space-Based Observatory 

Proposals currently exist for all-sky, space-based 
cosmic-ray detectors such as the Extreme Universe Space 
Observatory (EUSO) [TU] and the Orbiting Wide-Angle 
Lens (OWL) [H]. In addition, work is currently under¬ 
way to combine datasets from two large ground based 
experiments, the Pierre Auger Observatory (Auger) in 
the southern hemisphere, and Telescope Array (TA) in 
the northern hemisphere [12]. This paper will use EUSO 
as the example for a full-sky observatory, but our conclu¬ 
sions will apply to any full-sky observatory. 

EUSO is a down looking telescope optimized for near 
ultraviolet fluorescence produced by extended air showers 
in the atmosphere of the Earth. EUSO was originally 
proposed for the International Space Station (ISS), where 
it would collect up to 1000 cosmic ray (CR) events at and 
above 55 EeV (lEeV = 10^® eV) over a 5 year lifetime, 
far surpassing the reach of any ground based project. 

It must be emphasized that because previous data were 
so sparse at energies which would be accessible to EUSO, 
upper limits on anisotropy were necessarily restricted to 
energies below the threshold of EUSO. EUSO expects 
many more events at ~ 10^° eV, allowing an enhanced 
anisotropy reach. In addition, EUSO would observe more 
events with a higher rigidity TZ — E/Z, events less bent 
by magnetic fields; this may be helpful in identifying 
point sources on the sky. 

C. Space-Based Advantages 

EUSO brings two new, major advantages to the search 
for the origins of extreme-energy (EE) GRs. One advan¬ 
tage is the large held of view (FOV), attainable only with 
a space-based observatory. With a 60° opening angle for 
the telescope, the down pointing (“nadir”) FOV is 

7r(/iiss tan(30°))2 « ~ 150, 000 km^ . (1) 

We will compare the ability to detect large scale 
anisotropies at a space-based, full-sky experiment with 
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that of a ground based partial-sky experiment. For ref¬ 
erence, we will use the largest ground based cosmic ray 
observatory, the Pierre Auger Observatory [T3]. Auger 
has a FOV of 3,000 km^. Thus the proposed EUSO FOV, 
given in Eq. is 50 times larger for instantaneous mea¬ 
surements (e.g., for observing transient sources). Mul¬ 
tiplying the proposed EUSO FOV by an expected 18% 
duty cycle, yields a time averaged nine-fold increase in 
acceptance for the EUSO design compared to Auger, at 
energies where the EUSO efficiency has peaked (at and 
above ^ 50 — 100 EeV). Tilting the telescope turns the 
circular EOV given in Eq. into a larger elliptical FOV. 
The price paid for “tilt mode” is an increase in the thresh¬ 
old energy of the experiment. 

The second advantage of a space-based experiment 
over a ground based one is the coverage of the full-sky 
(47r steradians) with nearly constant exposure and con¬ 
sistent systematic errors on the energy and angle reso¬ 
lution, again attainable only with a space-based obser¬ 
vatory. This paper compares full-sky studies of possible 
anisotropies to partial-sky studies. The reach benefits 
from the 47r sky coverage, but also from the increased 
statistics resulting from the greater FOV. 

In addition to the two advantages of space-based ob¬ 
servation just listed, a third feature provided by a space- 
based mission may turn out to be significant. It is the 
increased acceptance for Earth skimming neutrinos when 
the skimming chord transits ocean rather than land. On 
this latter topic, just one study has been published El- 
The study concludes that an order of magnitude larger 
acceptance results for Earth skimming events transiting 
ocean compared to transiting land. Most ground based 
observatories will not realize this benefit, since they can¬ 
not view ocean chords, although those surrounded by ice 
such as IceTop m may also benefit from the ice as well. 

The outline of this paper is the following: We present 
the difference between partial-sky exposure and that of 
full-sky exposure in section [iTj In section |III| we review 
spherical harmonics as applied to a full-sky search for 
anisotropies, the power spectrum, and anisotropy mea¬ 
sures. In section [Tv| we explain the particular interest in 
the dipole {i = 1) and quadrupole {£ = 2) regions of the 
spherical harmonic space as well as the techniques used 
here and in the literature for reconstructing the first two 
spherical harmonics. We also discuss the difficulties in 
differentiating dipoles and quadrupoles with partial-sky 
coverage; these difficulties are not present in full-sky cov¬ 
erage. In section]^ we present the results of our analysis 
for a pure dipole or quadrupole. Finally, some conclu¬ 
sions are presented in section lYg 


II. COMPARISON OF FULL-SKY PROPOSED 
EUSO TO PARTIAL-SKY AUGER 

The Pierre Auger Observatory is an excellent, largest- 
in-its-class, ground-based experiment. However, in the 
natural progression of science, it is expected that even- 



FIG. 1. Auger’s exposure function normalized to / uj{Q)dfl = 
dvr. Note that the exposure is exactly zero for declinations 45° 
and above. 


tually ground-based observation will be superseded by 
space-based observatories. EUSO is proposed to be the 
first-of-its-class, space-based observatory, building upon 
ground-based successes. 

The two main advantages of a space-based observa¬ 
tory over a ground-based observatory are the greater 
FOV leading to a greater exposure at EE and the full- 
sky nature of the orbiting, space-based observatory. We 
briefly explore the advantage of the enhanced exposure 
first. The 231-event sample published by Auger over 9.25 
years of recording cosmic rays at and above ~ 52 EeV [16] 
allows us to estimate the flux at these energies. The an¬ 
nual rate of such events at Auger is ~ 231/9.25 = 25. 
For simplicity, we consider a 250 event sample for Auger, 
as might be collected over a full decade. 

Including the suppressed efficiency of EUSO down to 
~ 55 EeV reduces the factor of 9 relative to Auger down 
to a factor ~6 for energies at and above 55EeV. We arrive 
at the 450 event sample as the EUSO expectation at and 
above 55 EeV after three years running in nadir mode 
(or, as is under discussion, in tilt mode with an increased 
aperture but reduced PDM count). A 750 event sample 
is then expected for five years of EUSO running in a 
combination of nadir and tilt mode. Einally, the event 
rate at an energy measured by High Resolution Fly’s Eye 
(HiRes) is known to significantly exceed that of Auger. 
This leads to a five year event rate at EUSO of about 1000 
events. Thus, we consider the motivated data samples of 
250, 450, 750, and 1000 events in the simulations that 
follow. 

Now we turn to the 47r advantage. Auger’s exposure 
only covers part of the sky, and is highly nonuniform 
across even the part that it does see, as shown in Fig. 
The relative exposure is given explicitly as m 
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uj{6) cx cos ao cos 6 sin am + oim sin ag sin 6 
fo for^>l 

am = < TT for^ < -1 ^2) 

lcos“^^ else, 

cos 0m — sin ao sin S 

where f = ---, 

cos a cos 0 

and where a;((5) is the relative exposure at declination S, 
ao = —35.2° is Auger’s latitude, 0m = 80° is the new [T^ 
maximum zenith angle Auger accepts. We have assumed 
that the detector is effectively uniform in RA and that 
any variation (due to weather, down time, tilt of the ma¬ 
chine, etc.) does not significantly affect the uniformity of 
the exposure in RA. Auger recently modified their accep¬ 
tance from 0m = 60° —?► 80° with the extension calculated 
using a different metric: The S') 1000) technique is used 
for zenith angles 0 € [0°,60°], and the A^ig muon based 
technique is used for the new range, 0 € [60°, 80°]. These 
inclined events extend Auger’s reach up to a declination 
on the sky of -1-45°, as can be seen in Fig. In contrast, 
a space-based observatory such as EUSO would see in 
all directions with nearly uniform exposure. Of course, 
Auger is an existing observatory, while EUSO is but a 
proposal. 


III. TOOLS FOR ANISOTROPY SEARCHES 

A. Spherical Harmonics on the Sky 

As emphasized by Sommers over a dozen years ago HZ], 
a full-sky survey offers a rigorous expansion in spherical 
harmonics, of the normalized spatial event distribution 
where fi denotes the solid angle parameterized by 
the pair of latitude (0) and longitude (</)) angles, 

“ Jdn Nin) " £ II - (3) 

J ^ ' £=0 |m|<Z 

i.e., the set {Yfm} is complete. Y(f2) is the number of 
events seen in the solid angle U. The spherical harmonic 
coefficients, aim, then contain all the information about 
the distribution of events. The set {F^m} is also orthonor¬ 
mal, obeying 

J dn Yi,m, m Yl,m, m = 6l,l, 5m,m, ■ (4) 

We are interested in the real-valued, orthonormal Wm’s, 
defined as 

{ P^{x){y/2cos{m(j))) m > 0 

Pe{x) m = 0 , 

Pm{x){V2sm{m(j))) m<0 

(5) 

where, P^ is the associated Legendre polynomial. Pi = 
Pm=o i® i'ii® regular Legendre polynomial, x = cos0, 


and the normalization-factor squared is (f, m) = 

4-77 (£+m)l 

The lowest multipole is the i = 0 monopole, equal to 
the average full-sky flux and is fixed by normalization. 
The higher multipoles (^ > 1) and their amplitudes aim 
correspond to anisotropies. Guaranteed by the orthogo¬ 
nality of the Yim’s, the higher multipoles when integrated 
over the whole sky equate to zero. 

A nonzero m corresponds to 2 jmj longitudinal “slices” 
(jmj nodal meridians). There are ^ -I- 1 — jmj latitudinal 
“zones” (£ — |to| nodal latitudes). In Fig. we show the 
partitioning described by some low multipole moments. 
Useful visualizations of spherical harmonics can also be 
found in Ref. [TB]. The configurations with (£, — jmj) are 
related to those with (£,-|-|m|) by a longitudinal phase 
advance (/> => 0 -I- f, or cos 0 => sin (j). 


B. Power Spectrum 


The coefficients of the real-valued spherical harmonics, 
the aim's, are real and frame-dependent. To combat that 
problem, we use the power spectrum defined by 

( 6 ) 

m——i 


That the Ci should be rotationally (coordinate) invariant 
is not obvious. Certainly the spherical harmonics coef¬ 
ficients, the aim's given in Eq. are coordinate depen¬ 
dent. A simple rotation in the </> coordinate will change 
the sin (j), cos (/) part of the spherical harmonic for m ^ 0 
and a rotation in the 0 coordinate will change the asso¬ 
ciated Legendre polynomial part (P™(d)) for £ ^ 0. So 
only the £ = m = 0 monopole coefficient is coordinate 
independent. However, the power spectrum Ci is invari¬ 
ant under rotations. A recent derivation of this fact can 
be found in the appendices of Ref. (TH] . 

A simple approximation for the number of cosmic rays 
necessary to resolve power at a particular level is to count 
the number Nz{£,'rn) of nodal zones in each Wm- Each 
Yim has 


Nz{£,m) 


£+1 

2\m\[£ -I- 1 — jmj) 


m = 0 
m ^ 0 


(7) 


nodal zones. The average over m of the number nodal 
zones at a given £ is, 


(Nz(£)}= 3^^^^^\^ (2£^ + 4£ +3). (8) 


For low values of £, this returns the obvious results, 
(Nz(£ = 0)) = l,(Nz(£ = 1)) = 2. For large £, 
{Nz{£)) —t £'^/3. If we make the simple assumption of re¬ 
quiring 0{1) event per nodal zone to resolve a particular 
term in the power spectrum, then, for large £ we require 
~ £‘^/3 to resolve Ci. Thus, the rule of thumb is that our 
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FIG. 2. Nodal lines separating excess and deficit regions of sky for various (£, m) pairs. The top row shows the (0,0) monopole, 
and the partition of the sky into two dipoles, (1,0) and (1,1). The middle row shows the quadrupoles (2, 0), (2,1), and (2, 2). 
The bottom row shows the ^ = 3 partitions, (3, 0), (3,1), (3, 2), and (3, 3). 


EUSO fiducial samples of 450, 750, and 1000 events can 
resolve the C^’s up to an value of the mid-30’s, mid- 
40’s, and mid-50’s, respectively, i.e., (using 9 ~ can 
resolve structures on the sky down to 2-3°. A ground- 
based observatory, due to having fewer events and no 
full-sky coverage, would do much worse. We note that 
the statistical error in angle estimated here for EUSO is 
well-matched to the expected systematic angular resolu¬ 
tion error ~ 1° of EUSO. 


C. Anisotropy Measures 


Commonly, a major component of the anisotropy is 
dehned via a max/min directional asymmetry. 


a = 




G [0,1]. 


(9) 


A dipole (plus monopole) distribution is defined by a 
dipole axis and an intensity map given by /(O) oc 1 -I- 
Acos0, where 6 is the angle between the direction of 
observation, denoted by O, and the dipole axis. This form 
contains a linear combination of the Uim’s. In particular, 
a monopole term is required to keep the intensity map 
positive definite. One readily finds that the anisotropy 
due to a dipole is simply ud = A. 

A quadrupole distribution (with a monopole term but 
without a dipole term) is similarly defined, as /(O) oc 
1 — B cos^ 0. In the frame where the z axis is aligned 
with the quadrupole axis, the quadrupole contribution 
is composed of just the I 20 term. In any other frame, 
this I 20 is then related to all the Y 2 m’s, by the constraint 
of rotational invariance of the Cis mentioned above. In 
any frame, one finds that the anisotropy measure is ag = 
2 ^, and the inverse result is i? = 

Sample sky maps of dipole and quadrupole distribu¬ 
tions are shown in Fig. for both full-sky acceptance 
and for Auger’s acceptance, along with the actual and 
reconstructed symmetry axes. 


IV. RECONSTRUCTING SPATIAL MOMENTS 
A. Reconstructing a Dipole Moment 

Dipoles excite the specific spherical harmonics corre¬ 
sponding to Yim, with the three Yim’s proportional to 
X, y, and z. A dipolar distribution is theoretically mo¬ 
tivated by a single distant point source producing the 
majority of EECRs whose trajectories are subsequently 
smeared by galactic and extragalactic magnetic helds. 

With full-sky coverage it is straightforward to recon¬ 
struct the dipole moment so long as the exposure function 
is always nonzero (and possibly nonuniform). For the 
full-sky case (EUSO), we use the description described 
in m, which even allows for a nonuniform exposure, 
provided that the exposure covers the full-sky. 

Reconstructing any anisotropy, including the dipole, 
with partial-sky exposure is challenging. One approach 
for dipole reconstruction is that presented in [50]. We 
refer to this approach as the AP method. We note that 
this AP approach becomes very cumbersome for recon¬ 
structing the quadrupole and higher multipoles. Another 
approach for reconstructing any Yim with partial cover¬ 
age is presented in [21] , which we refer to as the AT-matrix 
approach. For the dipole with partial-sky coverage, we 
compare these two approaches to determine which opti¬ 
mally reconstructs a given dipole distribution. Our re¬ 
sult is seen in Fig. We consider 500 cosmic rays with a 
dipole distribution of magnitude OD^true = 1 oriented in a 
random direction. Using Auger’s exposure map, we then 
reconstruct the strength of the dipole using each method. 
This entire process was repeated 500 times. The recon¬ 
structed values oi au are a£),rec = 1.013 ± 0.101 and 
ttD.rec = 1.009 ± 0.084 for the AP and AT-matrix ap¬ 
proaches, respectively, where the uncertainty is statisti¬ 
cal. The mean angles between the actual dipole direction 
and the reconstructed dipole direction are 6 = 8.82° and 
6 = 6.41° for the AP and AT-matrix approaches respec- 
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FIG. 3. Shown are sample sky maps of 500 cosmic rays. The top row corresponds to the ao.true = 1 dipole, while the bottom 
row corresponds to the aqyrue = 1 quadrupole distribution. The left and right panels correspond to all-sky, space-based and 
partial-sky, ground-based coverage, respectively. The injected dipole or quadrupole axis is shown as a blue diamond, and the 
reconstructed direction is shown as a red star. We see that reconstruction of the multipole direction with an event number of 
500 is excellent for an all-sky observatory (left panels) and quite good for partial-sky Auger (right panels). In practice, an and 
aQ are likely much less than unity, and the event rate for EUSO is expected to be ~ 9 times that of Auger. Both effects on 
the comparison of Auger and EUSO are shown in subsequent figures. 


tively. The the results of the two approaches are compa¬ 
rable. Since the AT-matrix approach does slightly better 
than the AP method, we will use the AT-matrix approach 
for partial-sky dipole reconstructions in what follows. 

B. Reconstructing a Quadrupole Moment 

The physically motivated quadrupoles are the spherical 
harmonics corresponding to I 20 oc — T 20 represents 

an anisotropy that is maximal along the equator and min¬ 
imal along the poles (or, depending on the sign of 020 , the 
opposite). Such a distribution is motivated by the pres¬ 
ence of many sources distributed along a plane, such as is 
the case with the supergalactic plane. As a real, physical 
example of a well-known source-distribution distributed 
with a quadrupole contribution, we calculate the power 
spectrum as might be seen at the Earth for the 2MRS 
catalog of the closest 5310 galaxies above a minimum 
intrinsic brightness [52]. The catalog contains redshift 
information and contains all galaxies above a minimum 
intrinsic brightness out to z = 0.028 ~ 120 Mpc. As 
such, it is reasonable to suppose that EECRs come from 
these galaxies and, for simplicity, we implement uniform 
flux from each galaxy. 

In the left panel of Fig. we show the power spectrum 
that results for the known physical locations of these 
galaxies. In the right panel we show the power spectrum 


that results when each galaxy is weighted by the num¬ 
ber of events expected from it, i.e., by the inverse-square 
of the distance to the galaxy, a 1/d^ weighting. We re¬ 
mark that for the closest ~ 200 galaxies, the distance to 
each galaxy is known better from the direct “cosmic dis¬ 
tance ladder” approach than it is from the redshift, and 
we use these direct distance. For the farther galaxies, 
direct distances are less reliable, and we use the redshift- 
inferred distances. In this way, we also avoid any (possi¬ 
bly large) peculiar-velocity contributions to the redshifts 
of the nearer galaxies. 

It is instructive to compare the two panels. Without 
the 1/(P weighting (left panel), the intrinsic quadrupole 
nature of the distribution of 2MRS galaxies dominates 
the power spectrum; C 2 exceeds the other C^’s in the 
panel by a factor of > 5. In the right panel, galaxies are 
weighted by their apparent fluxes so the closest galaxies 
dominate. The large dipole is due to the proximity of 
Cen A, and the fact that the next closest galaxy, M87, is 
~ 4 times farther from the Earth. When determinations 
of the Cf’s are made, it is likely to be the dipole and 
quadrupole that will first emerge from the data based on 
the distributions of nearby galaxies. This quantifiably 
motivates our choice made in this paper to examine the 
dipole and quadrupole anisotropies. While the actual dis¬ 
tribution is likely a combination of dipole and quadrupole 
components, throughout this paper we consider the sim¬ 
pler cases where the distribution of sources has either a 
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FIG. 4. We simulated 500 cosmic rays with a dipole of amplitude of an,true = 1 pointing in a random direction, and with 
the Auger exposure. Then we reconstructed both the direction and the dipole amplitude 500 times with both reconstruction 
techniques (AP and A-matrix). In the left panel we show a histogram of the reconstructed values of ao.true, and in the right 
panel the angle (in degrees) between the correct dipole direction and the reconstructed direction. 


pure dipole anisotropy or a pure quadrupole anisotropy. 


As mentioned in section III C the quadrupole distri¬ 
bution that will be considered in this paper is of the 
form 1 — B cos^ 0, aligned with a particular quadrupole 
axis. The quadrupolar distribution is a linear combina¬ 
tion of the monopole term Iqo and I 20 oriented along 
the quadrupole axis. The distribution has two minima 
at opposite ends of the quadrupole axis and a maximum 
in the plane perpendicular to this axis. The quadrupolar 
data and the reconstruction of the quadrupole axis are 
shown in the lower panel of Fig. 

For the full-sky case, the method outlined by Som¬ 
mers in [m is used to reconstruct the quadrupole am¬ 
plitude and axis. It is possible to accurately reconstruct 
the quadrupole moment for experiments with partial-sky 
exposure at particular latitudes, independently of their 
exposure function. This is because there is very little 
quadrupole moment in the exposure function, as dis¬ 
cussed in Ref. m- By some chance, Auger is exactly 
at the optimal latitude in the southern hemisphere, and 
TA is very close to the optimal latitude in the north¬ 
ern hemisphere. Therefore we use Sommers’s technique 
for quadrupole reconstruction of both full-sky EUSO and 
partial-sky Auger. 


C. Distinguishing Between Dipoles and 
Quadrupoles 

One topic of concern is determining at what signifi¬ 
cance an injected dipole (quadrupole) distribution can 
be distinguished from a quadrupole (dipole), and from 
isotropy. Generally, the level of significance will de¬ 
pend on the number of observed cosmic-ray events, the 


strength of the anisotropy, etc. Fig. shows what hap¬ 
pens when Auger or EUSO attempt to reconstruct a pure 
dipole or a pure quadrupole when the signal is actually 
the opposite. The mean values and one standard devi¬ 
ation error-bars are derived from 500 repetitions of the 
given number of cosmic-ray events, where the dipole or 
quadrupole axis direction is randomly distributed on the 
sphere. The dashed lines in each plot are the 95% up¬ 
per limit for an isotropic distribution (i.e., atrue = 0). 
We see that as the actual anisotropy strength increases, 
quite a significant region of the parameter space would 
show an anisotropy in the absent multipole at the 95% 
confidence level when reconstructed by Auger. We also 
see that the relative size of the error bars reflects the 
statistical advantage of space-based observatories, while 
the central values of the data points, falsely rising with 
atrue for Auger but constant for EUSO, reveals the sys¬ 
tematic difference of partial-sky coverage versus full-sky 
coverage. 


This entire discussion is easily understood in the con¬ 
text of the “interference” of spherical harmonics which 
have been effectively truncated on the part of the sky 
where the exposure vanishes. The various truncated har¬ 
monics interfere heavily, a fact that is built into the K- 
matrix method (and into any method that attempts to 
reconstruct spherical harmonics based on only partial- 
sky exposure). Even though the true exposure of EUSO 
won’t be exactly uniform, the fact that it sees the entire 
sky with nearly comparable coverage means that the in¬ 
dividual spherical harmonics are non-interfering, and so 
can be treated independently. 
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FIG. 5. The power spectrum (see Eq. for nearby galaxies out to z = 0.028 (to d = 120 Mpc) based on their positions 
(left panel), and weighted according to 1/d^ (right panel). The 2MRS catalog [22] includes a cut on Milky Way latitudes 
|6| < 10° which is accounted for in the calculation of the power spectrum. C 2 is large because the galaxies roughly form a 
planar (quadrupolar) structure; Ci in the right panel is large because we are not in the center of the super cluster, thereby 
inducing a dipole contribution. (The relative scale between the ordinates of the two figures carries no information.) 


V. RESULTS 


In this section we tally our results. The standard pro¬ 
cedure involves simulating a number of cosmic rays with 
a given dipolar or quadrupolar anisotropy shape and am¬ 
plitude (atrue) aligned in a random direction. We then 
reconstruct the amplitude (arec) and direction (here we 
assume knowledge of the kind of anisotropy - dipole or 
quadrupole - unlike in section IV CI and compare to the 


true values. This process is repeated 500 times and the 
shown uncertainties are one standard deviation over the 
500 repetitions. 


A. Dipole results 

In Fig. we compare the capabilities of design EUSO 
and Auger to reconstruct a dipole anisotropy. In this 
comparison, both advantages of EUSO, namely the in¬ 
creased FOV and the Ltt sky coverage, are evident. The 
first panel shows how changing only the exposure func¬ 
tion between Auger and EUSO affects the value of the 
reconstructed dipole amplitude. For the same number 
of cosmic-ray events, the EUSO reconstruction is a bit 
closer to the expected value and has a smaller variation 
than does the Auger reconstruction. In the next panel 
we show the angular separation between the actual dipole 
direction and the reconstructed direction for Auger after 
a maximal amount of Auger data of ~ 250 cosmic-ray 
events, compared to EUSO’s minimal and maximal data 
sizes: 450 and 1000 cosmic-ray events, respectively. Even 
for a pure dipole. Auger will only reach 10° accuracy in 
dipole direction for a maximum strength dipole, ajj = 1, 
while EUSO does much better. In the third panel we 


compare both experiments at the same number of cosmic 
rays across a range of dipole strengths. Even if we assume 
that Auger will see significantly more cosmic rays than it 
is expected to, it still has a larger error in its ability to re¬ 
construct a dipole of any amplitude than EUSO. The low 
dipole magnitudes will always lead to a small erroneously 
reconstructed dipole due to random-walking away from 
zero. Finally, in the fourth panel we show the discovery 
power of each experiment to distinguish a dipole ampli¬ 
tude from isotropy. We see that Auger with 250 events 
would claim a discovery at five standard deviations above 
isotropy only if the dipole strength is 0.62 or greater - a 
situation that is unlikely given Auger’s anisotropy results 
to date [12]. EUSO could claim the same significance if 
the dipole amplitude is 0.37,0.30,0.27, or greater, for 
450, 750, or 1000 events, respectively. The EUSO signifi¬ 
cance should be enough to probe at high significance the 
weak signal currently reported by Auger. 


B. Quadrupole results 


In Fig. we again compare Auger and design EUSO 
in the context of quadrupole anisotropies. The same 
panels are plotted here as in Fig. except with an ini¬ 
tial quadrupole rather than dipole anisotropy, and a 
quadrupole reconstructed. We note that while the in¬ 
creased number of events that EUSO will detect will 
certainly lead to a better resolution of the quadrupole 
amplitude (as shown in the first and fourth panels) and 
direction (as shown in the second panel), we see that the 
full-sky coverage does not provide any benefit in this case 
(as deduced from the first and third panels). This result 
confirms the claims made in § IV B and in Ref. m- 
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FIG. 6. These panels show the results of attempting to reconstruct a dipole (quadrupole) when there is actually a quadrupole 
(dipole). The top two panels show the effect of attempting to infer a quadrupole moment from a pure dipole state of varying 
magnitudes while the bottom two panels show the effect of attempting to infer the dipole moment from a pure quadrupole 
state of varying magnitudes. The left two panels assume Auger’s partial coverage and 250 cosmic rays, while the right panels 
assume uniform exposure and the estimated number of events for EUSO (450 minimally, and 1000 maximally). The mean 
values and the one standard deviation error-bars are derived from 500 samplings. Note that the left most data point in each 
plot {ce{D,Q),true = 0) corresponds to the isotropic case, for which the dashed lines are the 95% upper limit. Finally, note that 
the vertical scales vary significantly between the partial-sky low statistics and full-sky larger statistics figures. 


Even though EUSO gains no benefit from its full- 
sky exposure for the determination of a quadrupole 
anisotropy, EUSO’s increased statistics will still lead to 
a detection sooner than Auger. Auger with 250 events 
would only be expected to claim a quadrupole discov¬ 
ery at five standard deviations above isotropy if the 
quadrupole strength is 0.67 or greater. EUSO could 
claim the same significance if the quadrupole amplitude 
is 0.47,0.36, 0.29, or greater, for 450, 750, or 1000 events, 
respectively. 


VI. CONCLUSION 

Many well motivated models predict, in the simplest 
limit, a dipolar or quadrupolar anisotropy in the EECR 
flux. The importance of the two lowest non-trivial orders 
(£ = 1,2) can be seen from the 2MRS distribution of the 
5310 nearest galaxies that was demonstrated in Fig. 
Due to the lack of any conclusive anisotropy from the 
partial-sky ground-based experiments, we explored the 
possible benefits that a full-sky space-based experiment, 
such as proposed EUSO, has over a ground-based exper¬ 
iment for detecting dipolar or quadrupolar anisotropies. 
In particular, we see that in addition to the increased 
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FIG. 7. Reconstruction of the dipole amplitude and direction across various parameters. Each data point is the mean value 
(and one standard deviation error-bar as applicable) determined from 500 independent simulations. The dipole amplitude and 
direction for Auger’s partial coverage were reconstructed with the A-matrix approach. The ordinate on the fourth panel, , 

labels the number of standard deviations above ao = 0. 


statistics that proposed EUSO brings over any ground 
based experiment, proposed EUSO significantly outper¬ 
forms Auger when reconstructing a dipole. Moreover, for 
inferences of both the dipole and the quadrupole, partial- 
sky experiments fail to differentiate between the two due 
to the mixing of the spherical harmonics when truncated 
by the exposure function. This situation is not present 
with all-sky observation, where the exposure function is 


nearly uniform and nonzero everywhere. 
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